Correlation energy of anisotropic quantum dots 
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We study the _D-dimensional high-density correlation energy of the singlet ground state of 
two electrons confined by a harmonic potential with Coulombic repulsion. We allow the harmonic 
potential to be anisotropic, and examine the behavior of Ec as a function of the anisotropy . 
In particular, we are interested in the limit where the anisotropy goes to infinity (a — >■ 0) and the 
electrons are restricted to a lower-dimensional space. We show that tuning the value of a from to 
1 allows a smooth dimensional interpolation and we demonstrate that the usual model, in which a 
quantum dot is treated as a two-dimensional system, is inappropriate. Finally, we provide a simple 
function which reproduces the behavior of E^ over the entire range of a. 
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I. INTRODUCTION 

The two-electron problem is one of the fundamental 
problems of quantum physics [THl] and, although it looks 
simple, it has only been solved in certain very special 
cases [51-fT5]. Many of the methods that have been de- 
veloped to provide approximate solutions to the two- 
electron problem have been central in the development 
of molecular physics and quantum chemistry |141 115j . 

The familiar Hartree-Fock (HF) model [16] treats a sys- 
tem as a separable collection of electrons, each moving 
in the mean field of the others. The HF solution pro- 
vides us with a good approximation to the energy and is 
widely applied to model complex molecular systems |17j . 
However, it is essential to understand its error [TB] 

Ec = E — i?HF, (1) 

which Wigner called the correlation energy [12]. Studies 
of correlation effects in two-electron systems are interest- 
ing in their own right, but also provide simple examples 
to test computational models [20] and shed light on more 
complicated systems [Hti^ . They have been extensively 
studied, for various confining external potentials, inter- 
acting potentials and degrees of freedom [51H27] . 

However, most previous studies have focussed on 
spherically symmetric external potentials, for anisotropy 
significantly complicates the mathematical analysis. 
This is unfortunate, for most real systems are not 
isotropic, and it is therefore important to understand how 
anisotropy affects the correlation energy. 



II. QUANTUM DOTS AT HIGH DENSITY 

Quantum dots are often modeled by electrons in a har- 
monic potential with Coulombic repulsion [55H3T]. Be- 
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cause experimental conditions strongly confine the elec- 
trons in one dimension, the model potentials are usu- 
ally spherical and two-dimensional. Calculations on such 
quantum dots have been used extensively in the devel- 
opment of exchange-correlation density functionals for 
low-dimensional systems in the framework of density- 
functional theory (DFT) [32H37]. 

In addition to experimental progress [55H1D] . many 
theoretical investigations have studied the effects of the 
confinement strength in the third dimension on the en- 
ergy of the quantum dot. These studies have used DFT 
[ill 132], HF [33], exact diagonalization [33] and exact 
solutions [I51 - ff7] . However, despite the importance of 
the correlation energy, only a few studies [48H50] have 
explored the confinement effect on E^- 

In this paper, we examine the effects of anisotropy on 
the energy of the nodeless ground state of two electrons in 
a D-dimensional harmonic potential, using atomic units 
throughout. We define the external potential by 
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where A governs the overall strength of the potential and 
aj e (0, 1] is the force constant for the Cartesian coor- 
dinate Xj. The isotropic case is obtained when all a^'s 
are equal. We are particularly interested in the behavior 
where one or more of the aj approach for, in such limits, 
the system is constrained towards a lower dimensionality. 

We restrict our attention to the high-density (A — > oo) 
limit [5Tti54] for it has been found that the high-density 
behavior of electrons is surprisingly similar to that at 
typical atomic and molecular electron densities [23H1Z]- 

The Hamiltonian describing this system is 

H = -^{\/l+\/l)+V{r^) + V{r2) + ^^, (3) 

where ri2 = [t-i — is the interelectronic distance and 
is the D-dimensional Laplace operator [55] |5B] . Seal- 
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ing all lengths by A, we obtain 



III. SECOND-ORDER ENERGIES 
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which is suitable for large-A perturbation theory with the 
zeroth-order Hamiltonian 
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and the perturbation 



ri2 



(5) 
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The Hamiltonian is separable and its eigenfunc- 
tions and energies are 



exp[-«,+4,)/(2a2)] 



X Hn^^^{xi.j/aj)Hn^^{x2,,/aj), (7) 



and 



(8) 



where = (nj^i, . . . ,ni^£i) holds the quantum numbers 
of the ith electron, Xij is the ith coordinate of the jth 
electron, and H„{x) is the nth Herniite polynomial |57j . 
Expanding the exact and HF energies as power series 



E = a2s(°) + Ai;(i) + + O(A-i), 



EnF = A^^F + + E^F + O(A-i), (10) 



one finds [52 EH [SI El] that 



= 4^), (11) 



and, therefore, that the limiting correlation energy is 



Ec = hm (E - Euf) = E'-^^ ~ E^ 



(2) 
HF- 
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In this paper, we show that i?c is strongly affected by 
the anisotropy and dimensionality of the potential. In 



Sec. Ill we use perturbation theory to obtain integral 



expressions for and 4f 

in an anisotropic quantum 



dot. In Sec. |IV[ we use the integral to express Ec as 
a infinite sum in a special case. Finally, in Sec. |Vj we 
present numerical results and discuss some of the impli- 
cations with regard to quantum dots and dimensional 
interpolation. 



The exact and HF second-order energies are [Ml ES] 
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Whereas the summation for the exact energy includes all 
states, the summation for the HF energy includes only 
singly-excited states 

Employing the Fourier representation 
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dk 
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one finds 
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(9) where F is the gamma function [5 7) . 



IV. 



CORRELATION ENERGY 



We now try to solve the integral 



(18) 



for special values of aj. The isotropic case, i.e. where all 
aj are equal, has been considered in detail in Ref. p4] . 
In the present paper, we generalize this to case where the 
aj take two distinct values. Without loss of generality, 
we let 5 be an integer such that < 6 < D, and set 



ai = ■■■ = as = a, 
as+i = ■ ■ ■ = ao = I. 



(19) 
(20) 



We assume a € (0, 1], i.e. the potential is stronger in the 
first S dimensions. 
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The integral (181 vanishes if any of the rii is odd. When 
all are even, it can be evaluated in hyperspherical coor- 
dinates EDI and one finds 
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where F is the Gauss hypergeomctric function |57j and 

D 

Since we only need to sum over terms where rij is even 
for all 1 < J < £), we can replace Uj 2 ?!, t hroughout. 
Substituting Eq. ^ into Eqs. ^ and ^ yields 

E,{a,D,S)- ^ 2 ^ 
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In the isotropic limit (a = 1), this becomes 
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as presented in Ref. [24] , 

In the anisotropic limit (a — > 0), we must consider 
two cases. If 5 = I? — 1, (23) becomes infinite, because 



the second-order energies and the correlation energy are 
unbounded for the one-dimensional dot [55 tl6TI - l63] . How- 
ever, if 5 < D — 1, we have 
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which confirms that, as the electrons are squeezed from a 
ZJ-dimensional space to a (D— 5)-dimensional space, their 
correlation energy tends smoothly toward the expected 
value for {D — 5) dimensions. 



TABLE I. Limiting correlation energies (miSh) for quantum 
dots with various a, D and 5. 
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FIG. 1. Variation of —Ec with D + (1 — ^/q)5, showing the 
correlation "pathways" that start from D = 3 (solid blue), 
D = 4 (dotted red), D = 5 (dot-dashed yellow), D = 6 
(dashed green) and D — 7 (solid black). 



V. NUMERICAL RESULTS AND DISCUSSION 

The correlation energies of quantum dots with various 
anisotropics a, for D = 3, ... ,7 and 6 — 1, . . . , D ~ 2 
are presented in Table |l] The energies for the spherically- 
symmetric states (a — 0, 1) coincide with the closed-form 
expressions reported in Ref. [57]. When 6 increases, the 
variation in correlation energy around a = 1/2 becomes 
flatter, but more pronounced near a = 0. In contrast, 
near a = 1, i?c remains independent of 6. Although not 
unique, pathways between dimensionalities are smooth. 
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FIG. 2. Variation of -E^{a, 3, 1) with a (solid). The fitting 
equation (261 is also shown (dashed). 



The correlation energy Ec{a, 3, 1) (in mE'h) can be ac- 
curately approximated using 



(a, 3, 1) « (co + cia + csa^) e"'^", (26) 



with 



Co 239.6, 
C2 = 557.5, 



ci = -360.3, 
C = 2.1736, (27) 



as shown in Fig. [2] By construction, this approximation 
is exact for a ^ and a = 1. 



VI. CONCLUSION 



and provide proper dimensional interpolations |55l 164) . 

Figure [T] shows the different correlation "pathways" 
when one starts from D = 3, . . . , 7 to reach 5 — 
1, — 2; the lower and the upper curves correspond 

to the smaller {5 — 1) and larger {5 — D — 2) values of 6, 
respectively. 

The D = 3 and 6 — 1 case models a quasi-2D quantum 
dot, confined in a three-dimensional harmonic well whose 
force constant is higher in one of the dimensions. How- 
ever, we see from Fig. [2] that, even when the anisotropy 
is large {a <^ 1), the correlation energy is significantly 
smaller than the D = 2 limit, which is frequently used 
to approximate the correlation energy of the quantum 
dot. Indeed, for a > 1/2, the correlation energy is much 
better approximated by the D = 3 limit. 

To illustrate the difference in correlation energy be- 
tween a strict and quasi 2D quantum dot, one can calcu- 
late Ec {a, 3, 1) for various values of a found in the liter- 
ature. In this way, we have found Ec = 73.8 for a — 0.43 
EHlllS], ^^c = 88.6 for a = 0.33 [Ml US] and E^ = 110.0 
viiE\, for a = 0.24 [39l|42]. Even for the smallest of these 
a values, the correlation energy is less than half of the 
2D limit (239.6 mEh). 



In this paper, we have studied the electron correlation 
of anisotropic quantum dots in the high-density limit. 
Using perturbation theory, we have solved the general 
Hamiltonian (|4]), and we have obtained integral expres- 
sions for i?*^^) and E^p. In the case where S dimensions 
are scaled by a factor of a with respect to the remaining 
(D — S) dimensions, we have expressed the correlation 
energy Ec as a infinite sum. 

Our numerical results reveal that Ec remains similar 
to that of the _D-dimension system for most a < 1, only 
increasing to that of the {D — (5)-dimensional state near 
a = 0. In the physically important D = 3 and 6 = 1 
case, the correlation energy is well approximated by the 
D = 2 limit only if a < 0.1. Such extreme anisotropy is 
probably difficult to realize experimentally. 
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